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Chapter  1:  Introduction 

Since  the  first  proposal  by  Kao  and  Hockham*  on  the  feasibility  of  optical  fiber  in  1966  and  the 
subsequent  demonstration  of  these  principles  at  Coming  Glass  Works  in  1970^  the  use  of  low 
loss  optical  fiber  has  rapidly  evolved.  The  interest  in,  and  massive  growth  of  fiber  optic  systems 
has  been  spurred  by  the  attractive  properties  offered  by  this  medium.  The  advantages  offered  by 
fiber  include  low  loss  transmission,  high  bandwidth,  small  size,  low  weight,  and  immunity  from 
external  electromagnetic  interference. 

Concomitant  with  the  growing  interest  and  deployment  of  fiber  optic  systems  is  the  development 
of  eomponents  to  be  used  in  these  systems.  As  a  result,  intense  interest  has  foeused  on  fiber 
Bragg  gratings  because  of  their  ability  to  be  used  in  many  different  applications.  These 
applications  include  rare-earth  doped  fiber  grating  lasers,^  dispersion  compensation,'*  wavelength 
division  multiplexing  (WDM),^  mode  eouplers,®  wavelength  stabilization  of  laser  diodes,’  hybrid 
fiber/semiconductor  lasers,*  fiber  amplifier  gain  control,®  grating  based  sensors,'®  time  delay  for 
phased  array  radar,"  and  nonlinear  effect  switches.*’ 

It  was  the  focus  of  this  work  unit  to  study  fiber  Bragg  gratings  and  establish  an  in-house 
capability  to  fabricate  fiber  gratings  for  RF  photonic  applications. 
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1.1  Historical  Background 


The  first  fiber  gratings  were  formed  external  to  the  core  region.  In  these  demonstrations  the  fiber 
cladding  was  side  polished  down  close  to  the  core  and  a  grating  was  introduced  by  either  placing 
a  metallic  grating  on  top  of  the  polished  region,'^  or  by  etching  a  grating  directly  into  the 
polished  surface.  These  gratings  interacted  weakly  with  the  evanescent  field  emerging  from  the 
core  region  and  were  challenging  to  fabricate.  As  a  result,  they  were  not  very  practical  for  real 
system  applications.  The  many  proposed  applications  inspired  by  these  gratings,  however, 
would  later  be  revisited  with  the  advent  of  in-core  fiber  gratings 

The  formation  of  a  permanent  grating  written  directly  into  the  core  of  a  fiber  was  first 
demonstrated  by  Ken  Hill  and  co-workers  at  the  Canadian  Communications  Research  Centre  in 
1978.'*  They  launched  high  intensity  488  nm  radiation  from  an  argon  laser  into  a  fiber  and 
observed  that  the  amount  of  light  back  reflected  out  of  the  fiber  increased  as  a  function  of  time, 
until  almost  all  of  the  light  transmitted  through  the  fiber  became  depleted  through  back 
reflection.  This  type  of  grating,  coined  a  “Hill  grating”,  was  formed  from  the  interference  of 
counter  propagating  beams  in  the  fiber  core.  Light  transmitted  through  the  fiber  experienced  a 
back  reflection  off  the  cleave  at  the  end  of  the  fiber.  This  back  reflected  light  interfered  with  the 
input  beam  and  the  resulting  standing  wave  pattern  was  recorded  in  the  fiber  core  through  a 
process  later  referred  to  as  photosensitivity. 
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The  formation  of  Hill  Gratings  was  an  unexplained  curiosity  imtil  the  important  findings  of  Lam 
and  Garside.*®  They  reasoned  that  a  two-photon  process  caused  the  grating  formation  and  that 
gratings  could  be  formed  more  efficiently  if  244  nm  radiation  was  used  instead.  This  inspired 
Meltz  and  co-workers  to  holographically  write  a  fiber  grating  with  244  nm  radiation  by  exposing 
the  fiber  core  to  two  interfering  beams  through  the  side  of  the  fiber.*’  By  varying  the  angle  at 
which  the  writing  beams  intersected,  fiber  gratings  could  be  made  to  act  as  Bragg  reflectors  at 
almost  any  wavelength. 

1.2  Basic  Concepts 

Since  the  initial  demonstration  of  fiber  gratings,  many  approaches  have  been  developed  to  write 
gratings  directly  into  the  core  of  the  fiber  (i.e.  holographic,  phase  mask,  prism,  etc.).  These  are 
all  similar  in  that  they  utilize  the  side  writing  technique.  They  vary  in  how  they  generate  the 
interference  pattern  used  to  form  the  grating. 

In  this  dissertation,  the  holographic  approach  was  used  because  gratings  of  any  fixed  period 
could  then  be  easily  created.  Fig.  1-1  illustrates  the  holographic  side  writing  technique.  A  single 
coherent  source  is  typically  split  into  two  beams  by  an  interferometer.  These  beams  are  made  to 
intersect  at  the  fiber  producing  a  periodic  interference  pattern  that  is  recorded  in  the  core  as  a 
refractive  index  modulation  with  period  A.  By  changing  the  angle,  frg,  at  which  they  intersect, 
the  period  of  the  interference  pattern  and  that  of  the  resulting  grating  may  be  varied. 
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Figure  1-1  Two  beam  holographic  exposure  of  the  fiber  core  using  the  side  writing  technique. 


A  signal  traveling  down  the  core  of  the  fiber  that  has  been  exposed  holographically  as  in  Fig.  1-1 
encounters  a  volume  hologram  where  the  writing  beams  overlapped.  This  signal  will  be 
diffracted  if  the  Bragg  condition  is  satisfied, 

2n^^Ksm.&  =  mXg 

where  n^^  is  the  effective  index  of  the  fiber  mode,  is  the  angle  of  incidence  on  the  grating,  A  is 
the  period  of  the  grating,  m  is  the  diffraction  order,  and  Xg  is  the  grating  wavelength.  Eqn.  1-1 
may  be  thought  of  as  a  phase  matching  condition  that  when  satisfied  implies  that  reflections  from 
successive  grating  planes  add  coherently. 
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Reflection  fiber  Bragg  gratings,  as  illustrated  in  Fig.  1-1,  have  fiinge  planes  which  are 
perpendicular  to  the  fiber  axis  (non-tilted  fringes).  These  gratings  act  as  a  wavelength  dependent 
filters  in  a  fiber.  When  the  wavelength  of  the  input  beam  matches  the  Bragg  condition,  the  angle 
of  incidence,  & ,  equals  the  angle  of  diffraction  and  the  beam  is  reflected  back  upon  itself  The 
grating  period  for  this  case  is  given  by, 

A  =  1-2 

Fiber  Bragg  gratings  can  be  characterized  by  their  transmission  or  reflection  spectra.  These 
spectra  may  be  measured  in  the  lab  with  the  aid  of  a  broadband  source  and  a  spectrum  analyzer 
or  calculated  by  use  of  the  coupled  mode  equations.  The  bandwidth  and  strength  of  these  spectra 
are  typically  proportional  to  the  length  of  the  grating  and  the  size  of  the  induced  index 
modulation. 

1.3  Organization 

This  report  is  composed  of  four  chapters,  including  this  introduction.  Chapter  Two  provides  the 
background  and  theoretical  models  used  to  predict  fiber  grating  performance.  In  Chapter  Three, 
a  discussion  of  fiber  grating  fabrication  is  presented  as  well  as  the  experimental  demonstration  of 
different  types  of  fiber  gratings.  Concluding  comments  are  made  in  Chapter  Four. 
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Chapter  2:  Theoretical  Performance  of  Fiber  Bragg  Gratings 

In  this  chapter  the  models  used  to  predict  the  theoretical  performance  of  fiber  Bragg  gratings 
are  presented.  Since  most  of  theory  is  based  on  mode  coupling,  Section  2.1  provides  a  brief 
introduction  into  the  basic  principles  of  optical  waveguides  used  throughout  this  ehapter. 
Section  2.2  deals  specifically  with  optical  fiber  waveguides  and  defines  the  modes  supported 
by  a  fiber.  Finally  in  Section  2.3,  these  coneepts  are  used  to  define  and  model  fiber  gratings. 

2.1  Dielectric  Optical  Waveguides 

An  optical  waveguide  is  a  structure  which  guides  light  from  one  point  to  another.  It  is 
typically  constructed  with  materials  of  different  dieleetric  constants  that  are  “stacked” 

f 

perpendicular  to  the  direetion  of  propagation  as  illustrated  in  Fig.  2-1 .  Light  propagation 
occurs  in  the  guiding  region  of  the  waveguide  through  total  internal  refleetion  at  the  material 
interfaces.  For  this  to  oeeur,  the  guiding  region  must  have  a  dielectrie  eonstant  (and  in  turn  an 
index  of  refraction)  of  greater  value  than  the  materials  surrounding  it. 

ric 

nf 

ns 

Figure  2-  1  Simple  slab  waveguide  with  eover,  film,  and  substrate  layers  and  nj>n^>n^. 

Light  propagates  in  the  longitudinal  direction  (z-direction)  in  the  film  layer  via  total  internal 

reflection  at  each  interface. 
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Following  Kolelnik’s  approach,'  it  is  assumed  that  the  electromagnetic  fields  in  a  waveguide 
have  a  periodic  time  dependence  and  can  be  written  in  the  form, 

^{x,y,z,t)  =  ^{x,y,z)e‘^  +  E\x,y,z)e~'‘^  2-1 

U(x,y,z,t)  =  W^x,y,z)e'‘"  +  ll\x,y,z)e~"^  2-2 

where  E(x,y,z)  and  H(x,y,z)  are  complex  amplitudes.  These  amplitudes  may  be  expressed  in 
terms  of  their  transverse  (x,y-direction)  and  longitudinal(z-direction)  components  relative  to 
the  waveguide, 

E(x,j;,z)  =  E,+E,  2-3 


H(x,y,z)  =  H,+H,. 


2-4 


Maxwell’s  equations  are  used  to  solve  for  the  electro-magnetic  fields  propagating  in  an  optical 
waveguide.  Assuming  a  lossless  media.  Maxwell’s  equations  can  be  written  as 

V  X  E(x,y,z)  =  -i(ajull(x,y,z)  2-5 

V  X  H(x,y,z)  =  ia)^(x,y,z) .  2-6 

In  solving  these  equations  boundary  conditions  must  be  satisfied  at  each  material  interface.  If 
we  assume  that  the  waveguide  in  Fig.  2-1  is  uniform  in  the  longitudinal  direction  and  changes 
only  in  the  transverse  direction,  the  resulting  solutions  must  satisfy  a  transverse  boundary- 
value  problem.  A  waveguide  mode  is  a  field  solution  to  this  problem  defined  as 


E(x,y,z)  = 

2-7 

H(x,y,z)=H,(x,y)e-'"^^ 

2-8 

where  the  subscript  v  denotes  the  mode  label  and  is  the  propagation  coefficient  of  the 
mode.  Every  distinct  mode  v  has  an  eigenvalue  and  a  unique  field  distribution. 
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In  general  there  are  two  basic  types  of  solutions,  guided  modes  and  radiation  modes.  Guided 
modes  are  a  set  of  solutions  that  are  confined  near  the  axis  of  the  waveguide  and  form  a  set  of 
discrete  /?’s.  Radiation  modes,  on  the  other  hand,  form  a  continuum  of  values  of  ft  and 
radiate  out  of  the  guiding  region. 

An  important  characteristic  of  the  modes  is  orthogonality.  This  means  that  modes  are 
independent  and  do  not  interact  with  each  other  as  they  propagate.  For  guided  modes  this 
requires[l]. 


where  P  is  the  power  carried  by  the  mode  and  5^^  is  the  Kroenecker  delta  function  (^^,  =1  for 
i=j:  ^^.=0  otherwise). 

An  important  consequence  of  mode  orthogonality  is  that  any  arbitrary  field  distribution  in  the 
waveguide  can  be  written  as  a  superposition  of  waveguide  modes.  For  an  arbitrary  field 
propagating  in  the  z-direction,  the  transverse  component  may  be  written  as 

2-10 

The  ability  to  represent  an  arbitrary  field  as  an  expansion  of  waveguide  modes  is  an  important 
tool  used  in  the  coupled  mode  formalism  of  Section  2-3. 
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2.2  Modes  of  an  Optical  Fiber 

An  optical  fiber  is  a  cylindrical  dielectric  waveguide.  It  consists  of  a  doped  silica  core 
surrounded  by  a  pure  silica  cladding.  Since  the  index  of  the  core  is  greater  than  the  index  of 
the  cladding,  guiding  is  achieved  via  total  internal  reflection  at  the  core/cladding  interface. 


To  calculate  the  modes  of  the  fiber,  the  wave  equation  is  solved  in  cylindrical  coordinates  to 
accommodate  the  geometry  of  the  fiber.  Due  to  the  complexity  of  the  solutions  for  this  case, 
it  is  easier  to  solve  for  the  longitudinal  field  components  first  and  then  use  Maxwell’s 
equations  to  find  the  remaining  transverse  components.  Following  the  treatment  given  by 
Pollock\  we  begin  with  the  scalar  wave  equation  for  the  longitudinal  field  component  E^, 


1  c/ 


2-11 


where  r  is  distance  along  the  fiber  radius,  ^  is  the  azimuthal  angle,  and  z  is  the  distance  along 
the  fiber  axis  (Fig.  2-2). 


Figure  2-2  Optical  fiber  waveguide  consisting  of  a  circular  core  surrounded  by 

a  lower  index  cladding. 
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is  defined  as  an  independent  fimction  of  r,  and  z,  respectively,  so  that  the  technique  of 
separation  of  variable  may  be  employed, 

E(r,  (l),z)  =  2-12 

where  the  e'"'  time  dependence  of  the  wave  has  been  omitted  for  clarity,  but  is  imderstood. 
Substituting  Eqn.  2-12  into  Eqn.  2-11  yields, 

+  4  +  kyR^e~‘^  =  0 . 

r  r 

Multiplying  Eqn.  2-13  by  rV R^e~'^  and  grouping  like  terms  results  in. 


2R  R 
r  — +r — 
R  R 


O" 

-r  a  +K„n  r  =  - — . 
^  "  O 


2-14 


The  left-hand  side  of  Eqn.  2-14  depends  only  on  R  while  the  right  side  depends  only  on  O . 
By  definition,  R  and  O  vary  independently.  Thus  both  sides  of  this  equation  can  be  set  equal 
to  a  separation  constant  v^: 


O"  +  v^O  =  0 , 


2-15 


and 


2  ^  R  oR  2/t22  o2  ^  \  n  n 
a-  dr  r 


2-16 


The  solution  to  the  differential  Eqn.  2-15  is  [2], 

[cos(v^ 
[sin(v^) 


rcos(v^) 


2-17 


Given  that  the  fiber  is  circularly  symmetric,  v  must  be  an  integer  to  satisfy  the  geometrical 
constraints  of  the  problem. 
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Eqn.  2-16  is  the  differential  equation  for  the  Bessel  functions.  Its  solutions  are  of  the  form/ 

2-18 


EXr,(p,z)  = 


+  c.c.  r  <  a 
\CKXyi')e''^ e~'^^  +  c.c.  r  >  a 


where  is  the  Bessel  function  of  the  first  kind  of  order  v,  and  is  the  Bessel  function  of  the 
second  kind  of  order  v.  k  is  the  transverse  wavevector  defined  as, 

K=-Jk^n  -0^  ,  2-19 

and  y  is  the  attenuation  coefficient  defined  as, 

2-20 

In  a  similar  fashion,  the  wave  equation  can  also  be  solved  for  resulting  in  [2], 


H,(r,(p,z)  = 


BJ^{Kr)e  +  c.c.  r  <  a 


2-21 


DK^{yr)e"''’’ +  c.c.  r  >  a 


The  transverse  field  components  can  now  be  determined  by  substituting  the  longitudinal  field 
solutions  into  Faraday’s  law  to  generate  E^,  ,  H,,  and  : 

2-22 


e-A' 


^SE,  (oudH^ 

,  B — ^+— — ^  I 

\  dr  r  dcp ) 


^  K 


h;=4 

K 


H  = 


-i 


\r  d(p  or  ) 


R—L 

K  d'  r  d(p) 

^pdH^  dE' 

— — -+(oe—^ 


\r  d(p 


d-  ) 


2-23 


2-24 


2-25 


12 


To  solve  for  the  guided  modes  in  the  core  (r<a),  and  from  Eqns.  2-18  and  2-21  are 


substituted  into  Eqns.  2-22  through  2-25.  This  results  in 
they  are  [2]; 

r  JXiai)  0  -K,(W) 


four  linear  equations.  In  matrix  form 


0 


2-26 


Jl 

aK^ 

icos^. 


JMa) 


K 


J,{m) 


ICO/J. 


Jl 

aK^ 


Jy  («a) 

Jy{m) 


pv 


ay 


2  ^yiW) 


-KM 

iCDJU 


r 


y  ay 


■K,  (^) 
KM 


A 

B 

C 

D 


=  0. 


For  non-trivial  solutions,  the  determinate  of  this  matrix  must  be  equal  to  zero.  Calculation  of 
the  determinant  yields  the  dispersion  relation,  or  characteristic  equation  of  the  fiber. 


L/  \  KJXica)  yKM 


Knl„J.  Q<a)  ^  Knl„,K,  {ya) 

kJM)  yKM 


2-27 


This  transcendental  equation  defines  the  allowed  eigenvalues,  B,  for  each  mode.  The  possible 
mode  solutions  for  a  fiber  are  labeled  as  HE^,  EH^n,,  TEq^,  or  Tmon,.  Each  mode  corresponds 
to  a  discrete  value  of  found  from  Eqn.  2-27.  The  mode  label  indices  v  and  m  are 
associated  with  the  transverse  variables  r  and  p ,  respectively,  and  are  thus  integers.  The  label 
V  is  defined  as  the  angular  mode  number  and  m  as  the  radial  mode  number. 


For  the  special  case  of  v=0,  the  eigenvalue  equation  can  be  simplified.  Only  TE  (Ez=0)  or 
TM  (Hz=0)  modes  exist  in  the  fiber  for  this  case.  Because  they  are  rotational  invariant 


(—  =  0 )  these  modes  propagate  as  meridional  rays  down  the  fiber.  When  v  does  not  equal 
a(p 


zero  the  solutions  to  Eqn.  2-27  are  the  hybrid  EH  or  HE  modes.  Hybrid  modes  consist  of  all 
six  field  components  and  thus  have  longitudinal  field  components  for  both  E  and  H.  The  EH 


13 


or  HE  designation  is  dependent  upon  which  longitudinal  component,  or  H^,  provides  the 
largest  contribution  to  the  mode.  Hybrid  modes  propagate  as  skew  rays  in  the  fiber  and  are 
rotationally  invariant.  In  terms  of  the  nomenclature,  the  lowest  order  mode,  or  fundamental 
mode,  of  a  single  mode  fiber  is  the  HEj,  mode.  This  is  followed  by  the  TEg,,  TMg,,  and  HEji 
modes.  As  discussed  below,  the  alternative  LP  mode  representation  is  often  used  in  place  of 
this  designation  because  it  is  easier  to  calculate  and  experimentally  observe. 

2.2.1  LP  Fiber  Modes 

Eqn.  2-27  can  be  difficult  to  evaluate,  especially  for  large  values  of  v.  If  it  is  assumed 
however  that  n^iaj  is  approximately  equal  to  n(,o„,  Eqn.  2-27  becomes  greatly  simplified.  This 
approach,  used  by  Gloge'*  in  his  weakly  guiding  approximation,  significantly  eases 
computation  but  still  yields  accurate  results  since  the  core/cladding  index  difference  of  most 
fibers  is  typically  very  small. 

Making  the  weakly  guiding  approximation  in  Eqn.  2-27  and  rearranging  terms  by  utilizing 
Bessel  function  identities  yields  [2], 

J,.^{Kd)  Ki_^{ya)  2-28 

Jj.,{Ka)  Kjiya) 

where  the  subscript  j  is  defined  in  terms  of  mode  number  v: 

j=l  TE,  TM  modes 

j=v-i-l  EHy  modes 

j=v-l  HEy  modes 
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The  basic  premise  behind  this  approximation  is  that  in  a  weakly  guiding  fiber,  linearly 
polarized  light  coupled  into  a  fiber  will  remain  linearly  polarized.  Polarized  modes  are 
obtained  mathematically  by  superimposing  certain  combinations  of  hybrid  modes.  Because 
the  resultant  superposition  modes  are  polarized  in  only  one  transverse  direction  (E,,  or  Ey), 
they  are  called  linear  polarized  or  LPj„  modes.  These  modes  are  degenerate  in  v  and  ^  and 
are  a  superposition  of  degenerate  EH^^.,  „  and  HE^.,  „  modes  and  degenerate  TE„„,  and  TM„„ 
modes. 

The  lowest  order  mode,  the  LPq,  mode,  corresponds  to  the  hybrid  HE,,  mode  and  is  the 
fundamental  mode  in  single  mode  fibers.  It  posses  a  two-fold  degeneracy  which  stems  from 
the  rkiad=ikore  approximation.  This  degeneracy  manifests  itself  as  two  orthogonally  polarized 
LPo,  modes  (E,,  or  Ey).  The  next  order  mode,  LP„  is  a  superposition  of  the  TEq,,  TMo„  and 
HE2,  modes.  This  and  all  other  higher  order  LP  modes  have  a  four-fold  degeneracy:  two 
modes  which  are  orthogonally  polarized  due  to  the  n<,,a<,=n<.„„  approximation  as  above,  and  two 
which  are  caused  by  a  dependency  in  the  solutions.  Calculated  intensity  plots  of  the  LPo, 
and  the  LP,,  are  provided  in  Fig.  2-3. 

In  general,  the  intensity  profile  of  a  mode  can  be  deduced  from  the  mode  subscripts  j  and  m. 
The  subscript  j  is  one  half  the  number  of  minima  (or  maxima)  in  the  intensity  distribution  as  a 
function  of  <p  from  0  to  2%.  Similarly,  the  subscript  m  is  the  number  of  maxima  in  the  pattern 
that  occurs  on  a  radial  line  from  zero  to  infinity.’ 
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Figure  2-  3  Calculated  intensity  plot  of  a  LPq,  fiber  mode  (top)  and  two  dimensional  cross 
section  of  a  calculated  LP,,  fiber  mode  (bottom). 
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2.3  Fiber  Grating  Caiculations 


The  basic  waveguide  and  fiber  concepts  presented  above  are  now  used  to  calculate  the 
spectral  performance  of  fiber  Bragg  gratings.  Central  to  the  approach  is  the  coupled  mode 
theory.  The  basic  coupled  mode  equations  are  derived  in  Section  2.3.1  using  the  formalism 
developed  by  Yariv.®  In  Section  2.3.2  these  results  are  applied  to  the  uniform  grating 
problem.  An  extended  coupled  mode  formalism  is  used  in  Section  2.3.3  for  the  case  of  non- 
imiform  gratings.  Finally  in  Section  2.3.4,  tilted  fiber  gratings  are  discussed. 

2.3.1  Coupled  Mode  Equations 

As  established  in  Section  2-1,  waveguide  modes  are  orthogonal  and  propagate  independently 
of  each  other.  In  an  ideal  waveguide  this  is  true  as  long  as  the  dielectric  constant  of  the 
guiding  region  remains  uniform  in  the  direction  of  propagation.  In  a  reality,  modes  can 
transfer  power  and  interact  with  each  other  due  to  deviations  from  ideal  conditions  in  the 
waveguide.  These  deviations  may  be  unintentional,  as  is  the  case  for  fabrication  errors  or 
bend  losses,  or  they  can  be  intentionally  caused,  as  is  the  case  for  photo-induced  fiber  Bragg 
gratings. 

To  calculate  mode  coupling,  deviations  in  the  waveguide  are  modeled  as  a  small  perturbation 
of  the  dielectric  constant  which  is  periodic  in  the  z-direction,  A£(x,y,  z) .  The  dielectric 
constant  is  then  written  as 

£ix,y,z)  =  {x,y)  +  A£{x,y,z) ,  2-29 

where  Sg{x,y)  is  the  ideal  or  unperturbed  part. 
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In  the  ideal  waveguide,  the  normal  modes  satisfy  the  wave  equation  for  the  unperturbed 
medium, 


(W^  +  e)E(x,y)=0, 


or  by  substituting  in  Eqn.  2-7, 

r 


+  ^  +  6)  jE„  (x,y)  =  0 . 

In  the  perturbed  waveguide,  on  the  other  hand,  As{x,y,  z)  can  couple  energy  from  a  given 
mode  into  other  modes.  As  discussed  in  Section  2-1,  the  resulting  field  in  the  perturbed 
waveguide  can  be  written  as  an  expansion  of  the  orthogonal  modes  of  the  unperturbed  system, 

E  =  (x,y)e"^"''  •  2-32 

m 

Substituting  Eqn.  2-32  and  Eqn.  2.29  into  the  wave  equation  (Eqn.  2-30)  yields, 

■^3,1  £X_,,»  1 

X  4-’  =0. 

[  +  co^ n{EXx,y)  +  Ae{x,y,z))A„, 

The  derivatives  of  E  with  respect  to  x  and  y  can  be  eliminated  by  use  of  Eqn.  2-3 1 .  Grouping 
the  remaining  derivatives  on  one  side  of  the  equation  and  rewriting  the  summations  in  terms 
of  integers  k  and  1  (which  can  only  assume  values  of  m)  yields. 


y  E^ix,y)e  ’'^^"  =-(0^nY,As{x,y,z)A^^{x,y)e'^'\ 

].  dz  dz  / 


Assuming  that  the  perturbation  is  small,  any  mode  coupling  that  occurs  will  be  slowly  varying 
and  thus  the  parabolic  approximation  can  be  made.  This  assumes  the  second  derivative  terms 
are  much  smaller  than  the  first  derivative  terms  and  may  therefore  be  neglected. 
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2-35 


S 

k 

Both  sides  of  Eqn.  2-35  may  be  multiplied  by  and  integrated  over  x  and  y, 


■i'^Pk^ 

az 


E*  (x,  ®  VS  (x,  y,  z)A,E,  {x,  y)e 


-iPiZ 


H2A^)JjEt(^,J^)-Ei(x,y)a4^= 

-oP/ii^A,  IJe*  {x, y)As(x,y, z)E, (x, 

This  allows  the  mode  orthogonality  relationship  to  be  used.  From  Eqn.  2-9  this  becomes 
^  J|E;(x,y)A£(x,y,z)E,(x,yy^^*'^'^<^4^.  ^  ^ 

Eqn.  2-37  can  be  written  in  terms  of  a  coupling  coefficient, 


lAirt-  ' 


where  is  defined  as 


2-38 


2-39 


The  coupling  coefficient  is  an  overlap  integral  whose  magnitude  describes  how  well  a  mode  k 
couples  to  a  mode  1  due  to  the  perturbation  As(x,y,z).  If  the  perturbation  As(x,y,z)  is  zero, 
the  coupling  coefficient  vanishes  and  no  mode  coupling  occurs.  If  As(x,y,z)  is  an  even 
function,  mode  coupling  can  only  occur  between  modes  of  the  same  parity  (i.e.  a  forward 
propagating  ZPqJ  mode  coupled  to  a  backward  propagating  LPQ■^  mode).  On  the  other  hand, 
when  A£(x,y,  z)  is  an  odd  function,  mode  coupling  can  only  occur  between  different  modes 
(i.e.  ->■  LP{[ ).  This  is  analogous  to  the  interaction  matrix  element  in  quantum 

mechanical  perturbation  theory.’ 
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2.3.2  Uniform  Fiber  Bragg  Grating 

A  uniform  grating  is  defined  as  having  a  sinusoidal  index  modulation  with  a  constant  or 
uniform  average  index  (Fig.  2-4).  For  this  model  it  will  be  assumed  that  the  uniform  fiber 
Bragg  grating  is  written  in  a  single  mode  fiber,  is  untilted,  and  is  a  reflection  grating.  In  this 
case  the  grating  will  couple  a  forward  propagating  mode  into  a  backward  propagating 
ZFo'j  mode.  A  grating  of  this  type  may  be  modeled  with  a  sinusoidal  index  perturbation  in  the 
fiber  core  given  by. 
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2-43 


'In 

M^)="co«+«lCOS(— z) 


where  A  is  the  period  of  the  grating  and  nj  is  the  magnitude  of  the  index. 


Figure  2-  4  Sinusoidal  index  profile  (Eqn.  2-43)  for  a  uniform  fiber  grating.  The  average 
index  is  uniform.  The  grating  has  a  length  L. 


The  first  step  in  modeling  a  grating  is  to  calculate  the  coupling  coefficient,  for  which 


bi.^x,y,z)  is  needed.  From  Eqn.  2-29, 


l^s{x,y,z)  =  E{^x,y,z)- s„{x,y). 


2-44 


where. 


s{x,y')  =  sj^{x,y'). 


2-45 


Applying  Eqn.  2-43, 


£(x,y,z) 


2  ,2n  .  2  2/2^  X 

z)  +  «,  cos  (— z) 
A  A 


2-46 
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where  n(,ore  is  the  index  of  the  fiber  core  and  n,  is  the  index  modulation  as  illustrated  in  Fig.  2- 

4.  Assuming  a  small  index  modulation  (ni«l),  the  last  term  of  Eqn.  2-46  will  be  very  small 

and  can  be  ignored.  A£(x,y,z)  is  then  approximated  by, 

A  /  ^  x2;r  ,  2-47 

A  s(x,  y,z)  =  2  cos(— 

A 

By  substituting  this  into  Eqn.  2-39  and  converting  to  cylindrical  coordinates,  the  coupling 
coefficient  becomes, 

0)  In 


where  the  subscripts  indicate  coupling  between  forward  and  backward  traveling  LPo,  modes. 
Notice  that  the  integral  on  the  right  side  of  the  equation  is  similar  in  form  to  the  Poynting 
equation.  Since  the  integration  is  performed  over  the  core,  it  represents  the  total  amount  of 
energy  propagating  in  the  core  of  the  fiber.  This  is  known  as  the  core  power  confinement 
factor,  r . 

The  coupling  coefficient  is  usually  written  in  reduced  form  as, 


<^ovor  =2x-cos(— z)  =  C*  , 

01+01 

where  k  is  defined  as  the  coupling  constant 


The  coupled  mode  equations  then  take  the  form. 


^  =  -iKA,{e 
oz 
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2-52 


dz 


=  iKA^{e  ^ 


-1P]2 


_/[2£+2^]z 

+  e  ^ 


) 


The  first  term  in  the  parentheses  of  both  equations  is  a  positive  exponential  that  is  rapidly 
oscillating,  and  hence  has  a  negligible  contribution.*  The  second  term  in  the  parentheses 
represents  the  phase  mismatch  between  the  two  modes.  The  grating  term  in  this  exponential 
can  compensate  for  the  mode  mismatch  when  the  Bragg  condition  is  satisfied.  Eliminating 
the  first  phase  term  yields, 

<^i  .  .  i2&  2-53 

— =  -ikA.b 
dz  ' 


dz  ' 


2-54 


where  is  the  Bragg  wavelength  of  the  grating,  X  is  the  wavelength  of  operation,  A  = 


2n 


is  the  grating  period,  and  S  is  defined  as  the  Bragg  detuning  parameter. 


c  ^  r,  2;r' , 

S  =  2\^~~  ~ 


1  1 


X  X 


■Bj 


2-55 


J  is  a  measure  of  the  Bragg,  or  phase,  mismatch  between  the  modes.  When  the  Bragg 
condition  is  met,  S  is  equal  to  zero  and  resonant  mode  coupling  occurs.  As  S  increases,  i.e. 
larger  Bragg  mismatching,  mode  coupling  decreases  to  zero. 


As  is  indicated  Fig.  2-4,  the  grating  is  defined  as  having  a  length  L  that  extends  from  z=0  to 
z=L.  If  we  assume  that  the  forward  propagating  wave  is  incident  at  the  start  of  the  grating  (let 
A,(0)=1)  and  there  is  no  energy  in  the  reflected  wave  at  the  end  of  the  grating  ( A2(L)=0),  then 
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a  closed  form  solution  subject  to  these  boundary  conditions  can  be  found  for  Eqns.  2-53  and 


2-54.  The  closed  form  solution  gives  the  reflected  power  R  = 


MO)' 

AiO) 
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„  (Klf  sinh^  i-JitcLf  -  (SLf )  2-56 

~  (fcLf  cosh^  -  (SLY)-  (SLf  ' 

Several  important  grating  properties  can  be  calculated  from  Eqn.  2-56.'®  For  example,  Fig.  2- 
5  illustrates  how  grating  reflectivity  changes  as  a  function  of  the  coupling  constant.  As  the 
fcL  product  increases,  the  reflectivity  increases,  the  grating  bandwidth  widens,  and  the 
strength  of  the  sidelobes  increases. 


Figure  2-  5  Reflection  spectra  of  a  uniform  fiber  grating  as  a  function  of  kL 
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The  bandwidth  of  the  grating  can  be  determined  by  manipulating  the  numerator  of  Eqn.  2-56. 
Dividing  the  equation  by  (kL)  ,  and  using  the  identity, 

sinh(/x)  =  zsin(x),  2-57 

the  numerator  of  Eqn.  2-56  becomes, 

sm{4{5Lf  -{Klf  .  2-58 

The  bandwidth  of  the  grating  is  defined  as  the  distance  between  the  first  zeros  in  the  reflection 
spectrum  on  either  side  of  the  peak  reflectivity.  These  zeros  occur  when  the  argument  of  Eqn. 
2-58  equals  ;r  or, 

AA  An  I  f  An  ^  ^  2-59 

Ag  ”o  1  ^AnlJ 

where  the  coupling  constant  from  Eqn.  2-50  has  been  used  Avith  An  =  Tn, .  There  are  two 

A 

important  limits  to  Eqn.  2-59.  The  first  occurs  when  An  is  small  (An  «  — ).  The  grating  is 

L 

considered  “weak”  and  in  this  limit  Eqn.  2-59  approaches, 

^^zern  _  2-60 

Thus  in  the  weak  grating  case,  the  bandwidth  is  inversely  proportional  to  length.  This  is  a 
reasonable  result  since  light  -will  traverse  more  grating  periods  in  a  longer  grating  and  have  a 
narrower  bandwidth  reflection  spectrum  as  a  result  of  phase  matching  over  a  greater  number 
of  fringe  planes.  A  plot  of  grating  bandwidth  versus  coupling  constant  for  a  weak  grating  is 

given  in  Fig.  2-6.  For  this  calculation  An  =  10”^  is  fixed  and  the  length  L  of  the  grating  is 

A 

varied.  The  grating  has  a  center  Bragg  wavelength  of  Ag  =  1 550n»i  and  satisfies  An«  —  . 
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Figure  2-  6  Calculation  of  grating  bandwidth  as  a  function  of  coupling  constant  for  the  case 
of  a  weak  grating.  The  length  L  of  the  grating  is  varied  for  several  cases  of  fixed  A« . 


For  the  opposite  limit,  A«  »-f,  the  grating  is  eonsidered  “strong”  and  Eqn.  2.59 

L 

approaehes, 

^  2-61 

In  this  limit,  the  bandwidth  becomes  independent  of  grating  length  and  depends  only  on  the 
index  modulation.  Light  propagating  through  a  strong  grating  is  mostly  reflected  before  it 
reaches  the  end  of  the  grating,  thus  negating  the  full  effect  of  the  grating  length. 
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The  bandwidth  of  a  strong  uniform  grating  as  a  function  of  the  coupling  constant  is  given  in 
Fig.  2-7.  An  is  varied  for  several  cases  of  fixed  length.  Again  =  I550nm  and  in  this  case 


An  » 


Figure  2-  7  Calculation  of  grating  bandwidth  as  a  function  of  coupling  constant  (An  is 
varied)  for  the  case  of  a  strong  grating.  Several  cases  of  fixed  length,  L,  are  plotted. 


Another  important  characteristic  of  a  grating  is  its  peak  reflectivity.  In  general,  maximum 
coupling  occurs  when  ^=0  (or  when  ^=1^).  When  resonance  does  occur,  Eqn.  2-56  reduces 
to 


^ax  =tanh^('^)- 


2-62 


In  Fig.  2-8,  is  plotted  as  a  function  of  the  coupling  constant.  As  shown  in  the  figure 
saturation  begins  to  occur  for  xi  >  4 . 


Figure  2-  8  Calculation  of  maximum  grating  reflectivity  as  a  function  of  coupling  constant. 

Eqn.  2-62  may  be  used  to  calculate  A«  from  experimental  data  if  the  grating  length  and 
reflectivity  are  knoAvn.  This  approach  is  used  in  Chapter  3  to  compare  these  calculations  with 
experimental  results. 

2.3.3  Non-uniform  Gratings 

The  intensity  profile  of  the  laser  beam  used  to  write  a  fiber  grating  determines  the  index 
profile  produced  in  the  fiber.  The  index  profile,  in  turn,  determines  the  spectral  characteristic 
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of  the  grating.  The  reflection  spectrum  of  a  grating  can  thus  be  tailored  by  apodizing  the 
writing  beams  to  create  a  non-uniform  grating.  As  opposed  to  imiform  gratings  which 
typically  have  a  reflection  spectra  with  strong  sidelobes  (as  was  shown  in  Fig.  2-5),  non- 
uniform  index  profile  gratings  can  be  created  with  the  spectral  characteristics  required  for  a 
specific  application. 


Non-imiform  gratings  may  be  modeled  by  a  raised  cosine  index  modulation  given  by 


11 


2n: 

n{z)  =  +  «i(^)(l  +  cos(— z)) 


2-63 


where  ni(z)  is  a  spatially  varying  envelope  function  of  the  induced  index  perturbation. 
Although  many  profiles  are  possible,  ni(z)  will  take  the  form  of  a  Gaussian  in  this  dissertation 
unless  otherwise  stated.  For  the  gratings  of  interest,  this  is  an  accurate  representation  of  the 
index  profile  since  the  gratings  are  formed  with  laser  beams  having  a  Gaussian  profile.  Thus 
n,(z)  is  given  by,'^ 

z'  ,  2-64 

n,  (z)  =  n,exp(-4  ln2 - 7) . 

FWHM^ 


Fig.  2-9  is  a  representation  of  a  Gaussian  grating  index  profile.  The  center  line  in  the  figure 
represents  the  spaced  average  index  across  the  grating.  As  a  grating  is  written,  this  index  will 
increase  until  saturation  occurs. 
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In  the  uniform  grating  model  of  Section  2.3.2,  the  average  index  remained  independent  of 
grating  strength  and  the  peak  of  the  grating  reflection  remained  centered  at  Xg .  In  the  non- 
uniform  grating  case  the  increase  in  the  average  index  with  exposure  causes  the  grating 
spectrum  to  red  shift  in  wavelength.  Since  the  index  change  is  larger  in  the  center  of  the 
grating,  this  portion  is  red-shifted  by  a  greater  amount  than  the  rest  of  the  grating.  As  a  result, 
it  reflects  light  at  a  longer  wavelength.  The  peak  of  the  grating  reflectivity  (associated  with 
the  peak  index)  is  thus  shifted  away  from  Xg .  The  weaker  parts  of  the  grating,  on  the  other 
hand,  are  not  red  shifted  as  much  as  the  center  of  the  grating  and  have  their  local  resonances  at 
shorter  wavelengths.  A  shorter  wavelength  is  transparent  to  the  center  of  the  grating  and  is 
not  reflected  by  it.  As  a  result,  effective  Fabry-Perot  cavities  with  resonance  peaks  on  the 
short  wavelength  side  of  the  peak  reflectivity  are  formed  [1 1]. 

This  situation  can  be  analyzed  using  the  extended  coupled  mode  formalism  of  Erdogen  and 
Sipe.’^  Bound-mode  to  bound-mode  coupling  in  non-uniform  gratings  with  either  straight  or 
tilted  fringes  may  be  calculated  with  this  method.  In  this  approach,  the  coupled  mode 
equations  become, 

^  =  iP[g /O-(z)  +  5  -  +  ipg^K{z)v{z)  ^ 

X  =  -W{gf<^(.z)  +  ^^]v(z)  -  iPg,K{z)u{z)  ^ 

dz  2p  az 

Where  u(z)  and  v(z)  are  the  forward  and  backwards  propagating  waves  in  the  fiber, 

<j(z)  represents  the  slowly  varying  Gaussian  background  index,  k{z)  is  the  amplitude  of  the 
index  modulation,  (f>  represents  the  variation  or  chirp  in  the  period  of  the  grating,  and  gf  and 
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gb  are  the  coupling  coefficients  for  forward  and  backward  scattering  from  the  fiber  grating, 
respectively. 


Unlike  the  simpler  case  of  the  uniform  grating,  a  closed  form  solution  to  Eqns.  2-65  and  2-66 
is  difficult.  Therefore  numerical  analysis  must  be  used.  A  program  was  written  in  MATLAB 
using  an  adaptive  step  size  Runge-Kutta  integration  algorithm  to  accomplish  this.  Runge- 
Kutta  is  typically  used  for  boundary  value  problems  with  initial  conditions  at  a  single  point. 
For  a  reflection  grating  the  initial  conditions  of  u(0)=l  (arbitrary)  and  v(L)=0  are  usually  used 
for  the  forward  traveling  and  reflected  waves.  Redefining  this  into  the  necessary  format  of  a 
single  point  boundary  value  problem  (as  opposed  to  a  two  point  problem),  the  initial  boundary 
conditions  of  u(L)=l  and  v(L)=0  were  used  and  integration  was  then  performed  backwards 


from  z=L  to  z=0.  The  grating  reflectivity  is  then  calculated  at  z=0  hy  R  = 


v(0) 


u(0) 


A  constraint  in  performing  this  calculation  was  that  the  Runge-Kutta  routine  in  MATLAB 
could  not  converge  with  equations  written  explicitly  in  terms  of  imaginary  numbers.  To 
overcome  this  problem,  the  equations  were  factored  into  real  and  imaginary  parts  and 
separated  via  the  substitutions, 

u->u  +iu  2-67 

V  ->  V  +  iv"  2-68 

where  a  single  prime  denotes  the  real  part  and  a  double  prime  is  the  imaginary  part.  This 
produces  four  equations  which  are  then  solved  simultaneously. 
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To  verify  this  approach,  the  simpler  uniform  grating  which  has  a  closed  form  solution  was 
used  as  a  test  case.  The  grating  spectrum  was  calculated  both  numerically  and  with  the  closed 
form  solution  of  Eqn.  2-56  and  the  results  were  compared.  This  was  repeated  for  a  variety  of 
index  modulations  and  grating  lengths  to  ensure  validity.  An  example  of  the  results  for  the 

case  of  a  1  mm  long  grating  with  an  index  modulation  of  approximately  A«  =  10"^  is  shown 
in  Fig.  2-10.  In  this  case,  as  well  as  the  others  tested,  the  results  of  the  numerical  analysis 
matched  those  derived  from  the  closed  form  solutions. 

With  the  numerical  technique  proven  to  be  accurate  for  the  uniform  case,  it  was  used  to 
calculate  the  solution  to  Eqns.  2-65  and  2-66.  The  results  for  the  case  of  a  straight  grating 
with  a  Gaussian  envelope  of  5  mm  (FWHM)  and  a  Bragg  wavelength  of  1550  nm  is  shown  in 
Fig.  2-11.  The  peak  of  the  grating  reflectivity  is  shifted  from  due  to  the  non-uniform  index 
distribution  as  described  above.  Also,  there  are  resonances  on  the  short  wavelength  side  of 
the  grating  peak  as  expected.  The  number  of  Fabry-Perot  resonances  are  a  function  of  grating 
strength  and  increase  with  increased  index  modulation. 

2.4  Summary 

The  tools  used  to  model  fiber  Bragg  grating  characteristics  were  presented.  Several  different 
grating  types  were  discussed.  The  techniques  used  to  fabricate  these  grating  types  are 
discussed  in  Chapters  3. 


Figure  2- 10  Reflection  spectrum  of  a  uniform  grating  calculated  by  both  the  closed  form  solution 
(solid  line)  and  numerically  (plus  marks). 
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Chapter  3:  Fiber  Bragg  Grating  Formation 

In  this  chapter,  the  fundamental  concepts  behind  the  experimental  fabrieation  of  fiber  Bragg 
gratings  are  presented.  To  describe  the  formation  of  fiber  gratings,  a  review  of  the 
photoinduced  meehanism  responsible  is  first  presented  in  Section  3.1.  This  photoinduced 
change  in  germania-silica  fiber  can  also  be  enhanced  with  hydrogen  gas.  Section  3.2 
discusses  the  important  control  parameters  for  hydrogen  loading  and  describes  a  loading 
system  designed  and  constructed  to  sensitize  fibers.  In  Section  3.3,  a  general  overview  of 
grating  writing  techniques  is  given.  The  formation  of  unslanted,  refleetion  fiber  and  the 
experimental  setup  used  to  fabrieate  them  is  presented  in  Section  3.4. 

3.1  Photosensitivity  in  Optical  Fibers 

Photosensitivity  in  optical  fiber  is  defined  as  the  ehange  in  the  refractive  index  within  the 
fiber  due  to  laser  exposure,  usually  in  the  UV.  Although  the  details  of  the  process  are  still  not 
well  understood,  one  popular  model  attributes  photosensitivity  to  germanium  related  defects 
in  the  fiber  core  as  deseribed  below. 

Optical  fiber  is  typically  made  of  glass  with  a  tetrahedral  structure  consisting  of  silicon  and 
oxygen  atoms.  In  this  matrix,  silicon  has  four  valence  electrons  and  oxygen  has  two.  Thus 
each  silicon  atom  is  singly  bonded  to  four  oxygen  atoms.  The  cladding  of  a  fiber  is  composed 
of  pure  silica  (SiOj)  glass.  The  core  of  the  fiber,  on  the  other  hand,  is  usually  doped  with 
germanium  to  increase  its  index  of  refraction  (this  allows  total  internal  reflection  to  occur  at 
the  core/eladding  interfaee).  Germanium,  like  silicon,  has  four  valence  electrons  and  can 
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replace  silicon  in  the  glass  structure.  This  substitution,  however,  produces  defects  when  the 
germanium  doesn’t  bond  properly.  The  resulting  defects  are  known  as  color  centers  and  are 
believed  to  be  the  origin  of  the  photosensitivity  effect. 


Figure  3- 1  Basic  GeO  defect  structures. 

Two  basic  kinds  of  germanium  oxygen  deficient  defects  have  been  proposed.  They  are 
illustrated  in  Fig.  3-1 .  The  first  is  shown  at  the  top  of  the  figure  and  is  known  as  a  neutral 
oxygen  vacancy  (NOV).'  A  NOV  occurs  when  a  germanium  atom  has  a  silicon  or  germanium 
atom  as  its  nearest  neighbor  instead  of  an  oxygen  atom.  The  second  defect  at  the  bottom  of 


the  figure  arises  when  a  germanium  atom  is  bonded  to  only  two  oxygen  atoms.  This  leaves 
the  germanium  atom  with  a  lone  pair  of  electrons  which  are  unbonded.^  An  important 
difference  between  these  defects  is  that  the  NOV  defect  can  be  bleached  by  irradiation  in  the 
240  nm  wavelength  band.  It  is  therefore  commonly  believed  that  the  NOV  defect  is  primarily 
responsible  for  photosensitivity,  but  the  role  of  each  type  of  defect  and  the  interaction  between 
them  is  still  not  well  understood. 


UV  spectra  of  3  mol%  gemumia  MCVD  core  before  and  after 
VV  exposure 

(i)  UV  absorption  of  MCVD  core  prior  to  UV  exposure 

(ii)  UV  spectrum  of  MCVD  core  after  JOmin  UV  exposure 

(iii)  UV  spectral  changes  induced  by  30min  UV  exposure 


Figure  3-  2  UV  absorption  spectrum  of  3  %  mole  germanium  doped  fiber  perform  before  and 
after  UV  exposure.  Taken  from  [3]. 
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Both  defects  do  exhibit  absorption  in  the  240  nm  wavelength  band.  Fig.  3-2  is  the  measured 
UV  absorption  spectrum  of  a  germanium  doped  fiber  perform  both  before  and  after  exposure 
to  UV  light  at  248  nm.^  The  key  point  of  the  figure  is  that  after  UV  exposure,  there  is  a 
decrease  in  the  absorption  in  the  240  nm  band.  This  bleaching  suggests  that  the  defect 
responsible  for  the  initial  absorption  has  been  eliminated. 

One  of  the  most  popular  models  used  to  explain  this  phenomenon  is  the  color  center  model.'* 
In  this  hypothesis,  it  is  believed  that  UV  light  ionizes  the  oxygen  deficient  color  centers 
causing  electrons  to  be  released.  The  freed  electrons  are  then  trapped  in  neighboring  sites, 
creating  new  types  of  color  centers  (i.e.  Ge(l)  and  Ge(2)).’  Hence  while  some  defects  are 
bleached,  others  are  created  and  the  net  result  is  a  local  redistribution  of  electrons  which 
modifies  the  UV  absorption  profile  of  the  glass  (Fig.  3-2).  A  change  in  absorption  can  then  be 
related  to  a  change  in  the  index  of  refraction  by  the  Kramers-Kronig  relation,® 

=  3-1 

71  •'  <y  -  0) 


where  the  P  stands  for  the  principal  part  of  the  integral  and  Aor  is  the  change  in  absorption. 
Thus  a  change  in  the  absorption  in  the  core  of  the  fiber  caused  by  exposure  to  an  UV 
interference  pattern  can  result  in  a  nonvolatile  modulation  of  the  index  of  refraction  in  the 
core  and  can  be  used  to  form  a  grating  in  the  fiber. 

The  magnitude  of  the  induced  refractive  index  change  varies  from  fiber  to  fiber  and  is  a 
function  of  germanium  doping  levels.  This  is  logical  since  it  is  the  germanium-related  defects 
which  are  believed  responsible  for  the  effect.  It  has  been  shown  that  absorption,  and  thus 
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index  change,  depends  linearly  on  germanium  concentration  up  until  the  onset  of  saturation.’ 
Saturation  is  both  a  function  of  energy  density  and  duration  of  exposure.  For  standard 
communication  fiber  which  is  doped  with  3-6  mole  %  GeOj,  the  UV  induced  index  change 
typically  does  not  exceed  10'“.* 

3.2  Hydrogen  Loading 

A  technique  developed  to  increase  the  photosensitivity  of  almost  any  germanosilicate  fiber  is 
“hydrogen  loading.”’  High  pressure  loading  of  hydrogen  gas  into  the  fiber  core  causes  nearly 
every  germania  site  in  the  glass  to  participate  in  the  formation  of  defects.  This  dramatically 
increases  the  magnitude  of  the  refractive  index  change  obtainable  and  allows  strong  gratings 
to  be  written  in  almost  any  standard  telecommunications  fiber.  Index  changes  as  high  as  10’’ 
have  been  reported.*® 

The  amount  of  hydrogen  which  may  be  introduced  into  a  fiber  is  a  function  of  temperature, 
time,  and  pressure.  The  influence  of  these  parameters  on  the  diffusion  rate  of  hydrogen  into 
the  core  of  a  fiber  can  be  modeled  by  using  a  Fickian  diffusion  equation.  This  models  a 
circular  cylinder  in  which  diffusion  is  everywhere  radial.  To  calculate  hydrogen  diffusion 
into  the  fiber  core,  the  fiber  is  treated  as  a  cylinder  with  a  radius  equal  to  the  fiber  cladding. 
The  normalized  concentration  of  hydrogen  at  the  core  of  a  fiber  with  radius  a  is  then  given 
by" 

C  2y  exp(-£>Q'’r)Jn(ra„)  3_2 

Q  cc„Jx{aa„) 
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where  D  is  the  diffusion  constant  of  hydrogen  in  silica,  t  is  time,  r  is  the  radial  distance  from 
the  center  of  the  core,  and  ra„  and  aa„  are  the  roots  of  the  Bessel  function  of  the  first  kind  of 
orders  zero  and  one,  respectively.  The  diffusion  constant,  D,  is  temperature  dependent  and 
has  been  found  empirically  to  be'^ 

D  =  2.83x10"''  exp(-4834/r)  /  sec .  •!_  o 


Figure  3-  3  Normalized  concentration  of  hydrogen  gas  diffused  into  the  fiber  as  a  function  of 
radial  position  in  the  fiber.  The  multiple  traces  represent  concentration  vs.  loading  time 
labeled  in  days. 
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Fig.  3-3  presents  a  series  of  plots  of  the  normalized  hydrogen  concentration  versus  radial 
position  in  the  fiber  at  a  temperature  of  294  K.  Multiple  plots  are  used  to  illustrate  how  the 
hydrogen  concentration  in  the  fiber  increases  with  time  (measured  in  day  increments.)  At  this 
temperature  it  takes  approximately  13  days  for  the  hydrogen  concentration  to  reach  95%  of  its 
equilibrium  value.  This  is  more  evident  in  Fig.  3-4  which  presents  the  data  in  the  form  of 
concentration  at  the  core  versus  time.  Increasing  the  temperature  at  which  the  hydrogen  is 
loaded  into  the  fiber  speeds  the  diffusion  process  and  decreases  the  loading  time.  For 
example,  elevating  the  temperature  to  300  K  reduces  the  loading  time  for  a  95%  normalized 
concentration  to  approximately  nine  days  as  shown  in  Fig.  3-4.  Although  increasing  the 
temperature  reduces  the  loading  time,  it  also  reduces  the  amount  of  hydrogen  molecules 
which  can  be  absorbed  into  the  fiber,  or  solubility  S.  Solubility  is  inversely  proportional  to 
temperature  and  is  given  by  [10], 

S  =  3.37exp(1041.66/  T)  ppm/atm.  2  3 

From  Eqn.  3-3,  increasing  the  temperature  results  in  a  lower  total  concentration  of  hydrogen 
in  the  fiber  at  a  given  pressure.  However,  a  higher  concentration  can  be  obtained  by  the 
loading  pressure.  The  equilibrium  concentration  of  hydrogen  at  the  core  of  the  fiber  is  equal 
to  the  solubility  times  the  loading  pressure, 

C„=S-Pppm.  34 

Although  greater  concentrations  can  be  obtained  with  increased  pressure,  the  mechanical 
tolerances  of  the  loading  system  impose  a  limit  on  the  maximum  pressure  which  can  be  used. 


43 


Normalized  Concentration  at  Core 


In  addition  to  the  obvious  concerns  of  using  a  high  pressure  gas  system,  hydrogen  gas  is 
explosive,  weakens  the  structure  of  metals  over  time,  and  can  be  volatile  under  pressure. 

Thus  the  optimal  balance  of  pressure,  temperature,  and  time  is  ultimately  dependent  upon  the 
system  used. 


Figure  3-  4  Normalized  concentration  of  hydrogen  in  the  core  of  the  fiber  as  a  function  of 
time  at  300  K  and  294  K. 
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3.2.1  Hydrogen  Loading  System 


A  hydrogen  loading  system  was  designed  and  constructed  to  sensitize  optical  fibers  as 
described  above.  The  system,  shown  schematically  in  Fig.  3-5,  consists  of  two  functional 
sections,  a  loading  assembly  to  house  the  fiber  and  a  gas  supply/delivery  system  to  pressurize 
the  loading  assembly.  The  entire  system  is  constructed  from  stainless  steel  pipes  pressure 
rated  for  7500  psi  and  high  pressure  Swagelok  fittings.  The  weakest  fitting  is  pressure  rated 
for  3000  psi  with  a  4: 1  safety  ratio.  The  system  is  typically  operated  at  2000  psi,  leaving  a 
1000  psi  safety  margin. 

The  loading  assembly  consists  of  a  loading  chamber  which  holds  the  fiber  during  loading  and 
a  water-based  heat  exchanger  used  to  regulate  the  loading  temperature.  The  loading  chamber 
is  simply  a  meter  long  stainless  steel  tube  0.25  inch  in  diameter.  Fiber  is  spooled  lengthwise 
on  a  rod  and  inserted  and  sealed  inside  this  tube.  The  heat  exchanger  is  formed  by  connecting 
a  one  inch  diameter  stainless  steel  tube  between  two  one  inch  T-fittings.  The  smaller  diameter 
loading  chamber  is  centered  inside  this  one  inch  wide  assembly  and  connected  to  it  at  the  ends 
of  each  tee  via  a  reducing  adapter.  This  creates  a  sealed  space  between  the  outside  of  the 
loading  chamber  and  the  inside  of  the  T-fittings  and  pipe.  The  bottom  of  the  T-fittings  are 
used  as  a  water  inlet  and  outlet  from  this  space  to  a  temperature  controlled  water  supply  as 
shown  in  Fig.  3-5.  Temperature  regulated  water  is  circulated  in  this  space  to  maintain  the 
temperature  in  the  loading  chamber  at  a  desired  value. 
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Figure  3-  5  Schematic  of  the  hydrogen  loading  system  consisting  of  two  functional  parts,  the 
gas  delivery  system  and  the  loading  assembly. 

The  gas  supply  system  uses  both  hydrogen  and  helium  source  tanks.  These  tanks  are 
connected  via  thermoplastic  hoses  (rated  for  5000  psi)  to  a  quick  disconnect  fitting.  The  quick 
disconnect  allows  the  loading  assembly  to  be  pressurized  and  then  disconnected,  isolating  the 
assembly  from  the  hydrogen  source. 
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Helium  is  used  first  to  safety  test  the  system.  Unlike  hydrogen  it  is  not  combustible,  but  like 
hydrogen  it  has  a  small  atomic  number  which  makes  it  a  valid  substitute  for  leak  testing.  The 
helium  gas  is  gradually  introduced  into  the  system  by  increasing  the  supply  pressure  in  steps, 
one  section  at  a  time,  using  various  safety  valves  in  the  system.  With  each  step  increase  in 
pressure  leak  tests  are  performed  using  both  a  high-pressure  leak  detector  fluid  and  a  helium 
“sniffer”  (a  mass-spectrometer  to  detect  trace  quantities  of  leaking  helium).  When  the  system 
is  fully  pressurized  and  verified  to  be  leak  free  and  safe,  it  is  de-pressurized  via  a  venting 
valve  and  then  re-pressurized  with  hydrogen  following  the  same  steps. 

The  system  is  typically  operated  at  2000  psi  to  allow  a  1000  psi  safety  margin  for  the  weakest 
fitting  in  the  system.  From  Eqn.  3-3,  the  equilibrium  solubility  of  hydrogen  in  the  fiber  is 
approximately  116  ppm/atm  at  a  temperature  of  294  K.  At  a  pressure  of  2000  psi  the  core  is 
loaded  with  1.6  %  hydrogen,  a  concentration  which  has  experimentally  proven  to  be  adequate 
for  writing  strong  gratings  in  most  fibers. 

3.3  Writing  Techniques 

There  are  several  different  approaches  to  writing  gratings  in  fiber.  Fig.  3-6  summarizes  the 
general  techniques. 

The  phase  mask  technique  is  illustrated  in  Fig.  3-6a.’^  A  phase  mask  is  a  silica  substrate  that 
has  a  grating  pattern  etched  onto  its  surface.  The  plus/minus  1  diffracted  orders  (the  zero 
order  is  typically  suppressed)  exiting  an  illuminated  phase  mask  interfere  to  form  a  grating  in 
a  fiber  placed  behind  the  mask.  The  advantages  of  the  phase  mask  technique  are  ease  of  use 
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for  the  mass  production  of  gratings,  easy  alignment,  reduced  stability  and  coherence 
requirements  compared  to  holographic  techniques,  and  the  ability  to  computer  generate 
arbitrary  interference  profiles  and  chirped  gratings.  The  main  draw  back  to  this  technique  is 
that  a  mask  can  produce  a  grating  at  only  a  single  specific  wavelength. 


Figure  3-  6  General  methods  used  to  create  fiber  gratings  utilizing  the  side  writing  technique. 
They  include  a)  phase  mask  technique,  b)  holographic  technique,  c)  prism  technique,  and  d) 
point  by  point  method. 


Thus  there  is  an  associated  cost  for  making  a  new  mask  each  time  a  grating  at  a  different 
wavelength  is  needed. 


Another  common  approach  to  writing  fiber  gratings  is  the  holographic  approach  illustrated  in 
Fig  3 -6b.  A  two  beam  interference  pattern  is  created  by  an  interferometer  and  recorded  in  the 
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fiber.  The  main  advantage  of  this  approach  is  that  the  period  of  the  interference  pattern  may 
be  tuned  to  produce  gratings  that  operate  over  a  wide  range  of  wavelengths.  Coherence  and 
stability  requirements  are  the  main  difficulties  encountered  with  this  method.  This  approach 
is  discussed  in  more  detail  in  Section  3.4. 

A  variation  of  the  holographic  approach  is  the  prism  technique  picture  in  Fig.  3-6c.  The  prism 
acts  as  a  common  path  interferometer  making  it  an  inherently  more  stable  arrangement.  The 
illuminating  laser  beam  is  incident  on  the  prism  hypotenuse.  This  input  is  oriented  so  that 
half  the  beam  is  reflected  off  the  back  face  of  the  prism  and  interferes  with  the  other  half  of 
the  beam  that  is  refracted  directly  to  the  bottom  of  the  prism.  In  essence,  the  grating  is  formed 
by  folding  a  beam  unto  itself.  One  drawback  to  this  approach  is  that  the  folded  beam 
produces  a  “half  Gaussian”  index  profile  in  the  fiber  which  can  create  strong  side  lobes  in  the 
reflection  spectrum. ''' 

The  final  technique  illustrated  in  Fig.  3-6d  is  the  point  by  point  writing  method.  A  writing 
beam  is  pulsed  as  the  fiber  is  translated  so  that  the  individual  grating  lines  may  be  written  into 
the  fiber.  The  input  beam  must  thus  be  tightly  focused.  Also,  the  process  is  slow  and  requires 
high  precision  translation. 

3.4  Reflection  Fiber  Gratings 

In  this  section,  the  fabrication  of  unslanted  reflection  fiber  Bragg  gratings  is  discussed. 
Gratings  are  formed  using  the  holographic  methods  as  discussed  below. 
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3.4.1  Grating  Period  Calculation 

In  general,  diffraction  from  a  periodic  structure  can  be  described  by  the  grating  equation, 

.  -  .  -  wA  3-5 

«sin^,  =«sinw  + - 

A 

where  is  the  angle  of  incidence  of  a  wave  on  the  grating,  6^  is  the  angle  of  the  diffracted 
wave,  A  is  the  grating  period,  and  m  is  the  order  of  diffraction.  The  grating  equation 
expresses  the  phase  matching  condition  that  results  from  the  boundary  conditions  imposed  by 
Maxwell  equations.  For  each  diffraction  order,  it  predicts  where  energy  will  be  scattered  and 
at  what  wavelength  constructive  interference  from  the  scattering  point  occurs.  For  the  special 
case  when  the  incidence  angle  equals  the  angle  of  diffraction,  the  Bragg  condition  is  satisfied 
and  resonant  coupling  of  energy  can  occur. 


The  grating  equation  may  be  applied  to  the  problem  of  diffraction  from  a  fiber  grating  by 
recognizing  that  the  mode  propagation  constant  p  is  equal  to  kne^  and  that  sin^ . 

Substituting  these  equalities  into  Eqn.  3-5  yields, 

,  ,  27m  3-6 

-I-  —  , 


or, 


3-7 


Assuming  a  non-slanted  reflection  fiber  grating  in  single  mode  fiber  that  couples  a  forward 
going  LPqi  mode  to  a  backward  propagating  LPq,  mode,  Eqn.  3-7  becomes. 
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which  is  the  Bragg  condition. 


Eqn.  3-6  may  be  expressed  pictorially  in  the  form  of  a  Bragg  Diagram.  The  diagram  provides 
a  convenient  approach  to  calculating  the  grating  period,  the  angle  of  diffraction,  and  the 
interbeam  angle  between  the  writing  beams  when  they  are  at  a  different  wavelength  than  the 
signal  beam.  Fig.  3-7  is  an  example  of  a  Bragg  diagram  for  the  case  of  an  unslanted  reflection 
grating.  In  the  figure,  K  represents  the  grating  vector,  is  the  wavevector  of  the  write 
beams,  and  is  the  wavevector  of  the  signal  beams.  Again  assuming  a  single  mode  fiber  and 
coupling  of  a  forward  going  LPqi  mode  to  a  backward  propagating  LPqi  mode,  the  magnitude 
of  the  k  vectors  are  both  equal.  By  simple  addition  of  the  k  vectors  in  the  diagram,  the  grating 
vector  is  found  to  be, 

K  =  2k^n ,  3-9 

or. 


To  form  a  grating  with  this  period  the  writing  wavelength  and  geometry  must  be  taken  into 
consideration.  The  desired  wavelength  of  operation  of  a  grating  is  typically  different  from  the 
wavelength  used  to  write  the  grating  and  thus  the  lengths  of  k^  and  Iq  will  be  different  (the 
most  commonly  used  propagation  wavelengths  are  1300  nm  and  1550  nm  while  the  grating  is 
written  at  244  nm.)  As  a  result,  the  interbeam  angle  formed  between  the  grating  vectors  of  the 
writing  beams  and  that  of  the  signal  beams  must  thus  be  different  if  Bragg  matching  is  to  be 


Figure  3-  7  Bragg  diagram  for  a  non-slanted 
reflection  grating. 

achieved.  Pictorially,  the  interbeam  angle  (in  the  fiber)  between  writing  beams  with 
wavelength  needed  to  achieve  Bragg  matching  at  a  signal  wavelength  is  equal  to  the 
angle  formed  by  and  the  perpendicular  bisector  of  K.  From  the  diagram, 


3.4.2  Experimental  Setup 


The  experimental  setup  used  for  fabricating  reflection  type  fiber  gratings  is  similar  to  a  basic 
Mach-Zehnder  interferometer  and  is  illustrated  in  Fig.  3-8.  The  244  nm  wavelength  laser 
energy  necessary  for  writing  gratings  in  fiber  is  produced  by  an  intracavity  frequency  doubled 


Figure  3-  8  Experimental  setup  used  for  producing  reflection  fiber  Bragg  gratings. 


Argon-ion  laser  (Coherent  model  FRED).  Light  from  the  laser  is  expanded  using  a  simple 
Galilean  telescope  and  is  directed  to  a  beam  splitter  using  a  series  of  UV  mirrors.  The  light  is 
then  split  along  two  paths  which  are  later  recombined  at  the  fiber. 

The  path  lengths  of  the  interferometer  were  equalized  within  the  coherence  length  of  the  laser. 
A  UV  grade  fused  silica  compensation  plate  was  inserted  in  one  arm  of  the  interferometer  to 
balance  the  mismatch  in  optical  path  lengths  between  the  two  arms  created  by  the  beam 
splitter.  Also,  an  additional  mirror  was  added  to  one  arm  of  the  interferometer  to  allow  the 
beams  in  both  arms  to  undergo  the  same  number  of  reflections  and  have  similar  portions  of 
the  two  wavefronts  interfere  with  each  other  in  the  fiber  plane. 

A  cylindrical  lens  is  used  to  focus  both  beams  along  the  axis  of  the  fiber  core.  This  produces 
a  greater  intensity  interference  pattern  on  the  core  and  reduces  exposure  time.  Typically  the 
cylindrical  lens  is  placed  before  the  beam  splitter.  This  eases  alignment  and  allows  both  arms 
of  the  interferometer  to  focus  the  beam  at  the  same  place.  This  approach  works  well  for 
hydrogen  loaded  fibers.  For  the  case  where  the  fiber  is  not  hydrogen  loaded  and/or  photo¬ 
sensitivity  is  decreased,  a  separate  cylindrical  lens  of  a  shorter  focal  length  is  used  in  each  arm 
of  the  interferometer.  A  shorter  focal  length  lens  positioned  directly  in  front  of  the  fiber 
produces  a  smaller  spot  size  on  the  fiber  and  results  in  greater  exposure  energy  at  the  core. 

The  two  turning  mirrors  in  the  interferometer  determine  the  Bragg  angle,  ,  and  the  period 
of  the  grating  (Eqn.  3-1 1).  Both  mirrors  are  mounted  on  precision  motorized  rotational  stages 
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with  .001°  angular  resolution.  These  stages  allow  precise  adjustment  of  the  Bragg  angle  for 
tuning  of  the  grating  wavelength.  To  obtain  the  proper  Bragg  angle,  a  reference  position  is 
first  set.  The  tuning  mirror  of  each  arm  is  positioned  to  reflect  the  light  coming  from  the 
beam  splitter  through  a  set  of  centered  irises  (Fig.  3-8).  The  light  is  thus  reflected  along  a 
path  that  is  parallel  to  the  normal  of  the  fiber.  From  this  position,  the  mirrors  are  then  moved 
by  an  angle  0^12  (since  the  reflection  angle  off  a  mirror  is  26 )  to  produce  a  beam 
intersection  angle  of  6^  with  the  fiber  normal  (as  shown  by  the  alternate  interior  angles  in 
Fig.  3-8.) 

The  focused  beams  intersecting  on  the  fiber  must  be  carefiilly  aligned  so  that  they  overlap  and 
focus  in  the  core  region.  Two  methods  are  used  to  aid  in  this  alignment.  The  first  is  to  use 
the  scattering  of  UV  light  off  the  fiber  to  help  determine  the  orientation  of  the  cylindrical  lens. 
Scattering  patterns  have  been  studied  in  depth  by  Presby  and  Marcuse  to  help  measure  fiber 
uniformity  and  index  during  fiber  manufacture.*^  Here  these  patterns  are  used  to  determine 
the  best  placement  of  the  cylindrical  lens  with  respect  to  the  fiber.  A  UV  sensitive  screen  is 
placed  behind  the  fiber  to  capture  both  writing  beams  after  they  diverge  from  the  fiber.  The 
fiber  produces  a  diffraction  pattern  and  a  shadow  of  itself  in  both  beams.  While  viewing  these 
patterns  on  the  screen,  the  height  of  the  cylindrical  lens  is  adjusted  until  the  shadow  is 
centered  in  the  diffraction  pattern.  Any  tilt/rotation  introduced  by  misalignment  of  the 
cylindrical  lens  is  then  adjusted  by  rotating  the  cylindrical  lens  until  the  sinc-like  diffraction 
patterns  are  parallel  to  the  fiber  and  the  shadow  region.  The  size  of  the  shadow  region 
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compared  with  the  entire  beam  profile  also  provides  a  rough  estimate  of  the  proper  placement 
of  the  fiber  with  respect  to  the  focal  plane  of  the  lens. 

The  second  alignment  technique  is  based  on  the  fluorescence  in  the  400  nm  region  the  fiber 
produces  when  irradiated  with  light  at  244  nm  [14].  As  the  light  is  focused  down  tighter  on 
the  fiber,  this  fluorescence  increases.  This  is  thus  used  to  position  the  lens  so  that  it  focuses 
the  light  in  the  fiber  plane.  The  blue  fluorescence  escaping  out  the  end  of  the  fiber  is 
measured  with  an  optical  power  meter  and  the  optimal  fiber  position  may  be  determined  when 
the  measurement  is  peaked.  These  two  alignment  techniques  may  thus  be  combined  to 
provide  course  and  fine  adjustment  of  the  cylindrical  lens  with  respect  to  the  fiber. 

Since  the  typical  period  of  a  reflection  fiber  grating  is  on  the  order  of  0.5  pm  or  less  (Eqn.  3- 
6),  system  stability  and  vibration  isolation  are  important  issues.  Hence  all  the  components  in 
the  interferometer  are  mounted  as  low  to  the  optical  table  as  possible  for  greater  fringe 
stability.  A  Plexiglas  case  also  covers  the  interferometer  to  isolate  the  setup  from  air  currents. 
As  a  check  on  stability,  a  Michelson  interferometer  using  a  HeNe  laser  was  constructed  within 
the  Plexiglas  box.  The  fringes  from  this  interferometer  are  projected  onto  a  screen  and  fringe 
movement  is  monitored  to  indicate  any  instability. 

A  fiber  grating  is  typically  measured  in-situ  as  it  is  being  written  using  a  broadband  source 
and  an  optical  spectrum  analyzer  (OSA).  The  spontaneous  emission  from  a  semiconductor 
optical  amplifier  (SOA)  is  used  as  the  broadband  source  in  these  experiments  (Fig.  3.9.)  Prior 
to  exposure,  light  from  the  SOA  is  coupled  into  the  fiber  and  then  back  out  into  the  OSA 
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(Ando  model  AQ-6315B  with  a  maximum  0.05  nm  resolution.)  With  this  configuration, 
changes  in  the  power  spectra  of  the  SOA  caused  by  the  formation  of  a  grating  are  observed  on 
the  OSA  during  exposure.  These  changes  indicate  the  wavelength  of  operation,  the 
bandwidth,  and  the  strength  of  the  grating  being  formed.  Gratings  of  any  desired  reflectivity 
can  be  produced  by  simply  terminating  the  exposure  when  the  grating  strength  reaches  the 
desired  value.  By  using  a  3  dB  coupler  as  shown  in  Fig.  3-7,  a  grating  can  be  viewed  either  in 
transmission  or  reflection  on  the  OSA. 


3.4.3  Experimental  Results 

In  this  section  experimental  grating  data  is  presented  for  unslanted  fiber  Bragg  gratings. 
Common  to  all  of  the  data  provided  is  the  measurement  technique.  The  data  is  recorded  using 
the  OSA  configured  to  measure  the  transmission  of  the  SOA  through  a  fiber  during  exposure. 
Fig.  3-9  is  a  reference  of  the  power  spectra  for  a  1500  nm  SOA  with  no  grating  present.  In  the 
data  below,  the  gratings  in  the  1550  nm  region  appear  as  a  notch  in  this  spectrum.  SOAs 
centered  at  1300  and  800  nm  were  also  used  to  write  gratings  in  these  respective  wavelength 
ranges.  Their  spectra  are  similar  in  shape  and  bandwidth  to  that  of  Fig.  3-9. 
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Figure  3-  9  Transmission  spectrum  of  a  1550  nm  SOA  used  as  a  broadband  source  to  measure 
fiber  gratings. 

3.4.3.1  Uniform  Grating 

As  discussed  in  Section  2.3.2,  a  uniform  grating  has  a  sinusoidal  index  modulation  with  a 
constant,  or  uniform,  average  index.  To  form  such  a  grating,  spatially  uniform  writing  beams 
are  required.  The  writing  beams  produced  by  the  laser,  however,  have  a  Gaussian  profile.  A 
uniform  grating  is  thus  created  by  first  expanding  the  laser  beams  and  then  placing  an 
adjustable  slit  in  each  arm  of  the  interferometer.  The  slit  is  used  to  limit  each  beam  to  a  0.8 
mm  diameter  spot  of  uniform  intensity.  Experimental  data  is  given  in  Fig.  3-10  for  a  uniform 
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grating  formed  in  this  way  in  hydrogen  loaded  Coming  SMF-28  single  mode  fiber.  Also 
included  in  the  figure  is  the  calculated  theoretical  spectrum  obtained  from  using  the  coupled 
mode  theory  of  Section  2.3.2. 


Figure  3- 10  Measured  and  calculated  transmission  spectra  for  Bragg  reflection  in  a  0.8  mm 
long  uniform  grating  with  kL  =  1 .577. 

To  perform  the  theoretical  calculation,  the  coupling  constant  of  the  grating  can  be  determined 
with  the  aid  of  Eqn.  2-62, 

^max  =tanh^('^)>  3-12 
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using  the  known  length  of  the  grating  and  experimentally  measured  peak  reflectivity. 

Once  K  is  determined,  it  may  be  substituted  into  Eqn.  2-56, 

^  (KLfsinh\^l(KLf-(SLf)  3-13 

(tcLf  cosh'(V(xi)^  -  (SLf )  -  (SLf  ' 

and  the  theoretical  spectrum  can  be  calculated  as  a  function  of  wavelength. 

Since  the  grating  is  uniform,  it  is  expected  that  the  peak  of  the  grating  will  coincide  with  the 
Bragg  wavelength  at  the  center  of  the  grating.  Indeed,  in  Fig.  3-10,  the  experimental  peak 
reflectivity  occurs  at  the  center  wavelength  1300.785  nm.  The  measured  peak  of 
approximately  8  dB  at  this  wavelength  yields  a  value  of  xi  =  1 .577 ,  or  A«  =  8. 1 36  x  1 0”'' 
using  Eqn.  3-12.  Substituting  this  value  of  tcL  in  Eqn.  3-13,  a  reasonable  fit  to  the  measured 
data  is  obtained.  As  another  check,  Eqn.  2.59  predicts  a  grating  bandwidth  of  .00 \2(Ag) for 
this  case  or  1 .56  nm.  This  is  very  close  to  the  approximate  1 .5  nm  bandwidth  of  the  measured 
data  as  shown  in  Fig.  3-10. 

3.4.3.2  Non-uniform  Grating 

Gratings  with  a  non-uniform  index  profile  described  in  Section  2.3.3  were  also  written.  These 
gratings  had  a  Gaussian  index  profile  formed  with  beams  approximately  4  mm  in  diameter. 

An  example  grating  formed  in  hydrogen  loaded  Coming  SMF-28  is  given  in  Fig.  3-11. 

Notice  that  as  expected,  this  grating  is  no  longer  symmetric  like  the  uniform  grating  in  Fig.  3- 
10.  The  grating  reflectivity  for  this  case  is  at  least  30  dB,  the  measurement  limited  by  the 


60 


signal  strength  of  the  SOA  and  the  noise  floor  of  the  OSA.  Also  evident  are  the  Fabry-Perot 
resonances  on  the  short  wavelength  side  of  the  grating  that  were  predicted  and  discussed  in 
Chapter  2.  Fig.  3-12  is  a  plot  of  the  calculated  spectrum  a  4  mm  long  Gaussian  profile  grating 

with  (j  =  7  X 10“^ .  The  plot  is  similar  in  structure  and  strength  to  the  experimental  data  and  to 
first  order  provides  a  reasonable  match. 


Figure  3-11  Measured  transmission  spectrum  of  a  4  mm  long  fiber  Bragg  reflection  with  a 
Gaussian  profile.  Fabry-Perot  resonances  are  present  on  the  short  wavelength  side  of  the 
grating.  The  measured  transmission  of  the  grating  is  limited  by  the  80  dB  noise  floor  of  the 
OSA. 
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Figure  3- 12  Calculated  transmission  spectrum  of  a  4  mm  FWHM  Gaussian  index  profile 
grating  with  cr  =  7xl0"^,  A  =  .004 ,  and  n(.,3d=l .447. 

It  was  also  predicted  in  Chapter  2  that  the  average  index  of  a  non-uniform  grating  increases 
with  exposure  shifting  the  center  Bragg  wavelength.  Fig.  3-13  illustrates  this  with  a  grating 
written  in  hydrogen  loaded  QPS  fiber  (QPS-PFBG-1355-T)  that  was  measured  at  two 
different  times  (4  minutes  apart)  during  exposure.  Not  only  does  the  eenter  wavelength  shift, 
but  as  the  index  modulation  increases  over  time,  the  grating  bandwidth  increases  as  predicted 
by  Eqn.  2-61  for  a  strong  grating.  The  grating  reflectivity  is  also  greatly  enhanced  with 
increased  index  as  expected. 
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Figure  3-  13  Shift  in  center  wavelength  of  fiber  grating  from 
827  nm  to  829  nm  with  exposure  time. 
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3.5  Summary 

The  basics  of  fiber  grating  formation  were  given  in  this  chapter.  This  included  a  discussion 
on  a  current  theory  used  to  explain  photosensitivity  as  well  as  the  effects  of  hydrogen  loading 
on  the  process.  The  experimental  apparatus  used  to  load  fibers  and  write  gratings  was 
presented.  Experimental  results  for  fiber  gratings  were  given  and  compared  to  theory. 
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Chapter  4:  Conclusion 


4.1  Summary 

It  has  been  the  foeus  of  this  projeet  to  study  fiber  Bragg  gratings  and  develop  an  in-house 
eapability  to  model  and  fabricate  them. 

Waveguide  and  fiber  concepts  were  discussed  at  the  start  of  this  report  to  provide  a 
foundation  for  understanding  how  fiber  based  grating  devices  function.  To  this  end, 
coupled  mode  theory  was  used  to  model  reflection  fiber  gratings.  The  effects  of 
experimental  parameters  on  grating  characteristics  were  modeled  for  both  uniform  and 
non-uniform  grating  profiles  using  this  approach. 

Fabrication  of  gratings  in  standard  communication  single  mode  fiber  required  an 
enhancement  of  the  fiber’s  photosensitivity  by  way  of  high-pressure  hydrogen  treatment. 
The  theory  behind  this  technique  as  well  as  the  experimental  system  used  to  treat  the 
fibers  was  discussed.  Once  treated,  the  fibers  were  fabricated  using  a  holographic 
approach.  The  performance  of  gratings  fabricated  by  this  method  was  then  compared  to 
theoretical  predictions  and  was  found  to  be  in  good  agreement. 

Fiber  grating  show  great  potential  for  use  in  RF  photonic  systems.  Other  in-house  efforts 
are  currently  examining  the  use  of  fiber  gratings  in  phased  array  radar  beam  steering 
systems. 
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